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The interaction between polyelectrolyte and ionic surfactant is of great importance in different 
areas of chemistry and biology. In this paper we present a theory of polyelectrolyte ionic-surfactant 
solutions. The new theory successfully explains the cooperative transition observed experimentally, 
in which the condensed counterions are replaced by ionic-surfactants. The transition is found to 
occur at surfactant densities much lower than those for a similar transition in non-ionic polymer- 
surfactant solutions. Possible application of DNA surfactant complex formation to polynucleotide 
delivery systems is also mentioned. 
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INTRODUCTION 



Solutions containing polyelectrolytes remain an out- 
standing challenge to physical-chemistry. Due to the 
long-ranged nature of the Coulomb force our understand- 
ing of this class of polymers is still quite rudimentary. 
This situation can be compared to the one that existed 
in electrochemistry of the turn of the century, before De- 
bye and Hiickel (DH) presented their, now famous, theory 
of strong electrolytes . The fundamental question that 
must be addressed by any successful theory of polyelec- 
trolytes concerns with the role played by the counterions. 
In this respect, the traditional theories of liquid state are 
not of great help, since most of the approaches based 
on resolution of integral equations come to a dead-end 
when the numerical schemes used to tackle these dif- 
ficult problems fail to converge. The scaling theories, 
which have been so successful in elucidating the proper- 
ties of nonionic polymers 0] have, so far, proven futile 
in the face of large number of length scales relevant for 
polyelectrolyte solutions. What seems to be lacking is a 
mean-field theory of polyelectrolytes similar to the one 
created by Debye and Hiickel for simple electrolytes, and 
Flory for non- ionic polymers. In our previous work we 
have attempted to construct such a mean-field theory for 
one special class of polymer solutions, the rigid polyelec- 
trolytes The constraint of rigidity allowed us to 
study the effects of electrostatic interactions decoupled 
from that of conformational structure of polyions. The 
theory has proven to be successful in elucidating various 
thermodynamic properties of rigid polyelectrolytes in the 
presence, or in the absence, of monovalent salt. In this 
letter we shall present a theory of rigid polyelectrolyte 
and ionic-surfactant solutions. 

The interaction between polymers and surfactants is 
of great practical importance in areas as diverse as col- 



loidal stabilization, polymer solubilization, mineral flota- 
tion and flocculation, as well as various aspects of molec- 
ular biology and biochemistry . In many practical ap- 
plications the polymers are dissolved in some sort of polar 
solvent, typically water, leading to monomer ionization. 
This situation is very common for biological systems. For 
example, in an aqueous solution, the phosphate groups 
of a DNA molecule become ionized, giving it a net nega- 
tive charge. Similarly the phospholipids, which compose 
the cell membrane, in aqueous environment acquire a net 
negative charge. The repulsion between the like-charged 
molecules makes the introduction of a polynucleotide se- 
quence into a cell a formidable challenge to molecular 
biologists. It has been observed, however, that in the 
case of binding by ionic surfactant dissolved in a polyelec- 
trolyte solution, the adsorption isotherms show a striking 
degree of cooperativity. This surprising phenomenon sug- 
gests that ionic surfactants or ionic lipids can be used as 
a "packaging" in order to deliver polynucleotides into liv- 
ing cells. Indeed, some recent experiments demonstrate 
that the cationic lipid reagents provide some of the best 
methods available for the gene delivery systems ||. 



THE MODEL 



The solution under consideration consists of anionic 
polyions, monovalent salt, and cationic surfactant, inside 
a volume V (see Fig.l). It is important to remember that 
the overall system is charge neutral, what implies that the 
negative charge of polyions and the positive charge of sur- 
factants is counterbalanced by an appropriate number of 
counterions (univalent cations), and coions (univalent an- 
ions), respectively. Furthermore, to simplify the analysis 
we shall assume that all of the counterions are identical, 
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whether they are derived from polyions or from disasso- 
ciation of monovalent salt. A similar approximation will 
be made in the case of coions. 

In order to study the interaction between an ionic sur- 
factant and a polyelectrolyte we resort to the simplest 
possible model. The rigid polyions, of density p p , are 
represented by cylinders of length L and diameter a p . 
Each polyion has a charge —Zq, uniformly distributed 
along the length of the cylinder. The spacing between 
each charged group is b = L/Z. The cationic surfac- 
tants, of density p s , are modeled as flexible chains of n s 
monomers each, with the head group carrying a charge 
+q. For simplicity we shall assume that each monomer 
is a sphere of diameter a c . The density of counterions 
(cations) is p COU nt = Zp p + p salt , while the density of 
coions (anions) is p CO ion — Psalt + Ps- Both the coions 
and the counterions will be modeled as hard spheres of 
diameter a c and charge ±q located at their centers. The 
solvent (water) will be represented by a uniform medium 
of dielectric constant D. 



P— Pcoion : (2) 

Pt = Ps - m B p P , (3) 

where p+ is the density of free counterions, p_ is the den- 
sity of coions, and pf is the density of free amphiphiles. 

THE THEORY 

The main task of the theory is to determine the charac- 
teristic cluster size, i.e., to find the values of ng and tub- 
In order to achieve this goal the appropriate Hclmholtz 
free energy has to be constructed and minimized. The 
free energy can be decomposed into an electrostatic and 
an entropic contributions. The electrostatic contribu- 
tion arises due to the polyion- counterion- surfactant, the 
polyion-polyion, and the counterion- coion- surfactant in- 
teractions. The entropic part is the result of mixing of 
various species ||f§. 
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FIG. 1. A polyion of (cylindric) diameter a v and length 
L 2> a surrounded by spherical counterions and coions of 
diameter a c , and flexible surfactant molecules. The charge 
spacing on the polyion is b = L/Z, and the radius of the 
exclusion cylinder is a = (a p + a c )/2. 

The strong electrostatic interaction between the 
polyions, the counterions, and the surfactants leads to 
formation of clusters each made of one polyion, ns coun- 
terions, and rriB surfactants. In what follows we shall 
neglect the effects of polydispersity in cluster sizes, since 
it can be shown not to significantly affect the final results 
HQ . The counterion and surfactant association with the 
polyions reduces the number of free entities. Charge con- 
servation implies 
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FIG. 2. The DNA-DoTAB binding isotherms, X = ~4k B T, 
Z = 440, and £ = 4.17 for DNA at room temperature. The di- 
ameters of the polyions and the counterions are 27 A and 7.04 
A, respectively. The size of surfactant molecule is n s = 13. 
The concentrations of DNA and of added salt are 2 x 10 _6 M 
and 18mM, respectively. Note that at the transition the con- 
densed counterions are replaced by the ionic surfactants. The 
transition is found to be of the first order with the vertical 
line locating the point at which two local minima of the free 
energy become equal. The solid circles are the experimental 
data from Ref. [13] . 

The polyion- counterion- surfactant and the counterion- 
coion- surfactant contributions can be obtained in the 
spirit of DH theory ||,|,|. As a counterion or a sur- 
factant associates with a polyion it neutralizes one of 
its charged groups. Hence, the effective charge per unit 
length of a cluster, made of ns bound counterions and 
mB bound surfactants, is a c i = —q{Z — ns — wb) / 'L . Let 
us fix one such cluster and ask what is the potential that 
it feels due to the electrostatic interactions with the other 
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entities. In order to answer this question it is necessary 
to solve the Poisson equation, V 2< i>^ c ^ = —Airp q /D. Due 
to the hard core exclusion, for r < a = [a p + a c )/2 the 
charge distribution can be approximated as, 

S(r) 



Pq 



(4) 



2tt r 

For r > a, in the spirit of DH theory, we shall assume 
that 

p q = -(Z -Tib - m B )qp p + qp + e~^ °' ^ + 

-qp-e +f, ** W W +qp+e- f> * W V , (5) 

where j3 = l/(k B T). Upon linearization, the Poisson- 
Boltzmann equation can be easily solved to yield 
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D D KaKi{Ka) 



4> 



r < a . 
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where (ko) 2 = Anpl/T* , pi = p + + p_ + p+, and the re- 
duced density and temperature are respectively p* = p;a 3 
and T* = Dak B T/q 2 , while K n (x) are the n-order mod- 
ified Bessel functions of second kind. It is important 
to recall that the linearization of the Poisson-Boltzmann 
equation is justified by the renormalization of polyion 
charge through formation of clusters Hf§|7]]i| . 

In terms of this potential, the electrostatic energy of a 
cluster is 
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(8) 



with A$( d ) = + (2a cl /D)\n(r/a), for r < a; 

A$( c/ ) = $^ t } , for r > a. That is, we subtract the 
logarithmic potential produced by a line of charge, since 
it will only contribute to the self energy of a cluster. 
The electrostatic free energy density, f = —F/V, for 
the polyion-counterion-surfactant interaction is obtained 
through the Debye charging process, where all the parti- 
cles are charged from to their final charge QJ^] . 
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The electrostatic correlational free energy arising from 
the interactions between the free counterions, coions, and 
free surfactants is obtained using the usual Debye-Huckel 
theory @,|, 
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For sufficiently large separations the effective electro- 
static potential of interaction between two clusters sepa- 
rated by a distance r is p],|i"o[ 
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exp (—ast) 
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(12) 



where is the angle between two complexes. The short- 
ranged nature of the effective cluster-cluster interaction 
allows us to write its contribution to the free energy as a 
second virial term, averaged over the relative angle sus- 
tained by two macromolecules, 
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The entropic (mixing) free energy is obtained using 
the ideas derived from the Flory theory [ jTl| . In general 
jent _ ff nt ) where ff nt is the entropic contribution 
of each specie i. For free counterions and coions, 



P± - P± In <t>± 



(14) 



where <f>± = (irp±/6)(a c /a) 3 are the volume fractions oc- 
cupied by free counterions and coions. For flexible sur- 
factant chains the entropic free energy is (ll]] 



Pft nt = pt-ptH^/n s } , 
where the volume fraction of surfactant is 
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Finally, for complexes made of one rigid polyion, n B 
counterions, and m B surfactants, we find 



(3fT{Pp) = P P -Pp^ 
with 



4> c i(Z + n B + m, B ) 



KP-n 



1 



A(a/L) 



(Z + n s m B + n B )Cc 



2 1 

+ -z(n s m B + n B ) ( — 



(17) 



(18) 



and, C c ;, the internal partition function of a {n B ,m B )- 
complex, 



(19) 



The trace is taken over all possible configurations of n B 
counterions and m B surfactants associated to a polyion. 
The occupation variables, a c (t) and <r s (t), are such that 
a c (t) = 1 if the monomer t of the polyion is occupied by 
a condensed counterion, and a c (t) = if no counterion is 
associated at t. The occupation variable, a s (t), behaves 
in the same way, but for an association with surfactants. 
The Hamiltonian can be written as 
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An implicit constraint is that each monomer can have 
either a counterion or a surfactant associated, but not 
both. We have also made a simplifying assumption that 
the only effect of counterion or surfactant association is 
a local renormalization of a monomer charge. Note that 
the first term of the Hamiltonian couples all the sites, 
since it is due to the long-ranged Coulomb potential. The 
second sum runs only over the nearest neighbors, and is 
related to the hydrophobic interaction of the hydrocar- 
bon tails. The configurations in which agglomerates of 
surfactant molecules form are energetically favored, i.e. 
the hydrophobicity parameter is negative, x < 0. 

Even this, seemingly simple one-dimensional sub- 
problem, is impossible to solve exactly due to the long- 
ranged nature of the Coulomb force. We will, there- 
fore, resort to a mean-field bound given by the Gibbs- 
Bogoliubov-Feynman inequality. Defining x = ub/Z, 
and y = tub jZ we find fl 
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where S = Z [ip(Z) — ?/>(l)] — Z +1, ip(x) is the digamma 
function, and £ = /3q 2 /(Db) is the Manning parameter 
|l2| . In our previous study we have numerically checked 
that this, indeed, is a good approximation 

Minimization of the total free energy, / = f ent -^fP cs _|_ 
jees _|_ jpp^ Tfjifa reS p ec t to and raj allows us to de- 
termine the characteristic number of bound counterions, 
n* B , and of bound surfactants, m* B . We shall compare 
the predictions of our theory with the experimental mea- 
surements on DNA dodecyltrimethylamonium bromide 
(DoTAB) system The DoTAB is a cationic sur- 

factant with an alkyl chain of twelve carbons. We can 
estimate the value of the hydrophobicity parameter, Xi 
as follows. Consider a micell or a monolayer composed 
of DoTABs. The hydrophobic energy required to take an 
alkyl chain of twelve carbons from bulk hydrocarbon to 
water is measured to be approximately 20ksT |15|. We 
can interpret this energy as derived from the favorable 
interaction between the adjacent surfactants. Since each 
surfactant inside a micell or a monolayer has five or six 
nearest neighbors, we can estimate x ~ — 4fcgT. Clearly 
this is only a rough estimate but it should be sufficient 
to explore the ramifications of the new theory. 
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FIG. 3. The DNA-DoTAB binding isotherms for various 
concentrations of added salt: (a) 5 mM, (b) 18 mM, (c) 40 
mM. The other experimental parameters are the same as in 
Fig. 2. Note the change in the order of transition as it passes 
from continuous to discontinuous with an increase in the den- 
sity of monovalent salt. 

In discussions of adsorption it is traditional to define 
binding fractions, (3 C = n* B /Z and f3 s = m* B /Z. In Fig. 2 
we present the binding isotherms of DNA with dode- 
cyltrimethylamonium bromide, and compare it with the 
experimental data of jyj). It is evident that the agree- 
ment is quite good, without any fitting parameters! We 
note, however, that at densities of monovalent salt used 
in experiment, our theory predicts a first-order transi- 
tion, while the experimental data is more consistent with 
a second-order transition. In Fig. 3 we demonstrate that 
as the concentration of monovalent salt is lowered the 
transition becomes continuous. 



CONCLUSION 

We have presented a mean-field theory of 
polyelectrolytc-ionic-surfactant solutions. Although 
quite simple, our theory manages to capture the essen- 
tial physics of the problem. The most nontrivial aspect 
of polyelectrolyte-ionic-surfactant complex formation is 
that it occurs at extremely low densities, about a factor 
of twenty lower than the critical micell concentration 
(CMC) of pure amphiphile. This should be compared to 
the interaction of nonionic polymer with surfactant, in 
which case the binding transition happens at densities 
close to the CMC. Our theory explains this dichotomy 
in terms of strong electrostatic interactions which, in ad- 
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dition to hydrophobic forces, govern the polyelectrolyte- 
surfactant complex formation. 
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